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1. INTRODUCTION 

A common approach to the many-electron problem (atoms, molecules and 
solids) is its transformation into a fewer-particle problem. In Density Functional 
Theory 1,2 (DFT) one rewrites the ground-state energy of a many-electron system 
as a functional of just the one-electron density n(r), the diagonal part of the one- 
particle reduced density matrix. In recent years, much attention has been devoted 
to approaches which rewrite the system energy as a functional of the one- or the 
two-particle reduced density matrix. 3-9 With respect to DFT, these approaches have 
the advantage that the kinetic energy functional is known and that they provide 
more information about the many-body wavefunction. In this contest, the relevance 
of two-electron wavefunctions (geminals) has been pointed out: the use of an an- 
tisymmetrized geminal power wavefunction, i.e., a many-body wavefunction built 
in terms of geminals with the correct symmetry under particle permutation, pro- 
vides the basis for a formally correct one-particle reduced density matrix functional 
theory. 3 

With these concepts in mind, it is of great interest to start with a simple many- 
body system, the uniform electron gas, a limit which should be recovered by any 
approximate approach to the many-electron problem of nonuniform density. A more 
detailed discussion about reduced-density matrix functionals applied to the uniform 
electron gas can be found in Refs. 7, 10-12. In this paper, an approximation for the 
unknown two-electron wavefunctions (geminals) of the uniform electron gas is found, 
starting from the effective screened Coulomb potential proposed by Overhauser. 13 
The short-range part (small electron-electron distances) of the corresponding pair 
density is found to be in accurate agreement with the newest Quantum Monte 
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Carlo data. 16 This means that the short-range part of these geminals is reliable and 
could serve as a useful term of comparison for other two-electron approaches to the 
uniform electron gas. 7, 10-12 



2. ELECTRON PAIR-DENSITIES, ETC. 

Given an TV-electron wavefunction \I/(ri(7i, ...tnctn), we define the pair density 

N „ N 



p 2 (r,r') = N(N-l)Y^ I f[dv j \^(ra 1 ,r'a 2 , ...v N a N )\ 2 , 

<7( = 1 J j = 3 



(1) 



the one-particle density matrix 

N 

Pi(r,r') = N / **(rai,r 2 <72, ...rjva-jv)*(rVi, r 2 a- 2 , ...r N <r N )dr 2 ...dr N , (2) 

<7i=i •* 

and the one-electron density 

n(r) = Pl (r,r) = / dr'p 2 (r,r'). (3) 

While n(r)dr is the probability of finding an electron in dr, p 2 (r, r')drdr' is the 
probability of finding one electron in dr and another in dr' . We also define the 
pair-distribution function g(r, r'): 

p 2 (r, r') = n(r)n(r').g(r, r'). (4) 

By integrating Eq. (1) over r', we find that 

J dr'n(r')[g(r,r')-l] = -l. (5) 

In other words, the density n(r')[g(r, r') — 1] of the exchange-correlation hole around 
an electron at r represents a deficit of one electron. 



3. UNIFORM ELECTRON GAS (JELLIUM) 

The three-dimensional jellium model consists of N electrons enclosed in a box 
of volume V (periodically repeated in space) in the presence of a neutralizing back- 
ground of uniform positive charge density n + = N/V. The non relativistic jellium 
is thus governed by the hamiltonian (in Hartree atomic units): 



N 2 1 N I 

^EiHEi^t + a, ( 8) 
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where A represents the effect of the background. The hamiltonian of Eq. (6) de- 
scribes a model solid whose positive ionic charges are smeared throughout the whole 
crystal volume to yield a shapeless, uniform positive background (whence the nick- 
name of jcllium). When studying this model, one is usually interested in its macro- 
scopic properties, i.e., the thermodynamic limit (N, V — > oo keeping n = N/V con- 
stant) of its extensive physical quantities per particle or per volume. Two numbers 
are enough to describe its zero-temperature phase diagram, namely, the one-electron 
density n — N/V and the spin polarization ( = \Nf — N±\/N, where A'fm is the 
number of spin- up (down) electrons (N — Nf + Ni). Instead of the particle density 
n, it is often convenient to use the Wigner-Seitz radius r s (in units of the Bohr 

radius) given by r s = (4^ n) ; it is also useful to define the Fermi wavevector 

k F , simply related to the r s parameter by k F = (37r 2 n) 1//3 = (^ L ) 1 ^ 3 ^-- 

Beeing the system homogeneous and isotropic, the pair-distribution function 
only depends on r = |ri — r 2 |, and paramctrically on r s and £. 



4. JELLIUM FROM TWO-ELECTRON WAVEFUNCTIONS 

The pair-distribution function of jcllium can be built starting from two-electron 
wavefunctions. 13, 14 We first rewrite the non-interacting gas (ideal Fermi gas) in 
terms of two-electron wavefunctions. This first step can seem redundant, since 
the ideal gas can be written in terms of one-electron wavefunctions, as a Slater 
determinant of plane waves. However, treating first the ideal gas is essential to set 
the proper normalization. The two-electron wavefunctions for the interacting gas, 
in fact, are found by solving a scattering problem in an effective potential, with the 
normalization condition that the ideal gas is recovered when the potential is set to 
zero. In this paper, only the ( = case is analysed. For the generalization to the 
£ 7^ gas see Rcf. 14. 

4.1. Ideal Fermi gas 

If we select a pair of electrons at random in the spin-unpolarized uniform gas, 
there is one chance in four that they will be in the singlet state, Tl — IT) an d three 
chances in four that they will be in one of the triplet states, TT, II, Tl + IT- in 
the case of no electron-electron interaction, the corresponding two-electron spatial 
wavefunctions can be rewritten in the center-of-mass reference system as 

*(r,R) = ^e iK R (e jk r ±e- lk r ) , (7) 

where "+"is for the singlet state and "— " is for the triplet state, and 

R=|(ri+r 2 ), r = r 2 -n, K = kx+k 2 , k=|(k 2 -ki). (8) 

Beeing the system isotropic, it is convenient to expand the plane waves into spherical 
harmonics 

oo 

e lk ' r = ^(2^+l^P,(cos#)j,(fcr), (9) 
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where Pi are Legendre polynomials and ji are spherical Bessel functions. Then 

oo 

* singlct (r,R) = V2e lKR £ (2£ + 1) i e P e (cos9) Je (kr) (10) 



e=o 
even t 



*tripiot(r,R) = V2e tKR J2(2£+l)i i P e (co S 9)j e (kr). (11) 



i=i 

odd I 



The pair-distribution function can be obtained by ^/singlet and ^triplet by giving 
them the proper weight, 

5M = |(|* s i„ gl ct(r)| 2 ) + |(|* tri piet(r)| 2 ), (12) 

and by considering that in the uniform electron gas there is a probability p(k) that 
two electrons have a given relative momentum k = ||ki — k 2 |. In Eq. (12) the 
symbol () means, in fact, that an average over p(k) and over the solid angle has to 
be performed. If one is interested in the spin-resolved pair-distribution functions, 
<?lf (r) and g-\i{r), corresponding to parallel- and antiparallel-spin interactions, and 
such that for the unpolarized gas 

.9 = |(.9TT + .9Tt): ( 13 ) 

one has to consider that | of the triplet state (fj, + J,f) contributes to the antiparallel- 
spin correlations and | of it (|| and H) to the parallel-spin correlations. So, we 
have 

g n (r) = i(|*si„ g iet(r)| 2 ) + i(|* tripl ct(r)| 2 ) (14) 
3TTW = (l*tripict(r)| 2 ), (15) 

where () denotes again average over p(k) and over the solid angle. Performing the 
spherical average over the solid angle, we obtain: 



g n (r) = £(2^+l)(j 2 (fcr)) (16) 

e=o 

oo 

5TT (r) = 2^(2^+l)0' 2 (fcr)). (17) 



i=i 

odd l 



Equation (16) immediately gives the exact result for a noninteracting gas, i.e., 
.9Tl( r ) = 1 f° r each r. To obtain the noninteracting g^(r) from Eq. (17), we need 
to average over k. In the noninteracting electron gas, the probability distribution 
p(k) for k = i|k 2 — ki| can be obtained geometrically by considering two three- 
dimensional vectors ki and k 2 with < | ki(2) | < kp, where kp is the Fermi 
wavevector. The probability p(k) is then 

p(fc) = 24^-36^ + 12^, (18) 

rup rvp rvp 

with k ranging from to fcp (see Fig. 1). By using two known series which involve 
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Figure 1: Relative-momentum (k = i|k 2 — ki|) probability distribution p(k) for the 
nonintcracting gas. 



3l 



we can rewrite g^f as 



e=o 



sin(2a;) 
2x 



(19) 



g n (r) = J2W+mUkr)) £(-l) / (2^ + l)<j?(*r)> 



sin(2fcr) 
2kr 



(20) 



This means: 



r-kp 

Jo 



1 - 



sin(2fcr) 
2fcr 



2 v - 3 V 



A- 5 

'-V 1 ' M 



Performing this integral gives the known Hartree-Fock (noninteracting) g^, 



(21) 



.gf T F = 1-9 



sin(fci?r) — kpr cos(kpr) 



{k F rf 



(22) 



Numerically, in the range < kpr < 6 and with a truncation of the infinite sum 
over I at £ max = 7, Eq. (17) reproduces the exact Eq. (22) within an accuracy 
of 1CP 6 . When the scaled variable kpr (or equivalently r/r s ) is used, the pair- 
distribution function of the noninteracting gas does not depend on r s . The explicit 
dependence on r s only appears when Coulomb repulsion is taken into account in 
the wavefunction. 



4.2. Interacting electron gas 
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The interacting case can now be treated by introducing an effective potential 
which describes the electron-electron interactions in a uniform electron gas. The 
spherical Bessel functions je which appear in Eqs. (16) and (17) are solution of the 
noninteracting radial Schrodinger equation 

ut = ui = krji(kr). (23) 

If we introduce an effective radial potential V(r, r s ) which depends on the electron- 
electron distance r and parametrically on the electron density r s , we can solve the 
corresponding interacting radial Schrodinger equation, 

ui = u e = kr R e (r, k, r s ). (24) 

The radial functions Re will depend parametrically on k and on r s . We can insert 
them into Eqs. (16) and (17), and find the corresponding and g^: 

oo 

9n(r,r s ) = £(2*+l)<i$(r,A;,r a )) (25) 

1=0 

oo 

g n (r,r s ) = 2 £ (2£ + l)(i? £ 2 (r, k, r s )). (26) 

odd I 

To compute the average over all the possible relative k (represented again by the 
symbol () ) one should in principle use the interacting momentum distribution, which 
deviates from the Fermi step function of the ideal gas. This would slightly change 
p(k) of Fig. 1 by adding a "tail" for k > k F and lowering the maximum. For prac- 
tical purposes, the use of the noninteracting p(k) is enough to give good results. 14 
Notice that unless the potential V(r,r s ) is very sophisticated, the treatment just 
described will fail to recover the long-range correlations, which are mainly governed 
by collective modes, and will fail to satisfy the particle-conservation sum rule of 
Eq. (5). 



dr 2 



1(1+1) 



+ fc 2 



*-^-V(,r.) + *. 



5. SOLUTION OF THE OVERHAUSER MODEL 

In this section, we compute an interacting pair-distribution function following 
the procedure just described by using the simple model potential V(r, r s ) proposed 
by Overhauser. 13 This simple model gives very good results for the short-range 
(r < r s ) part of g(r). 14 

5.1. The Overhauser potential 

Overhauser 13 proposed a simple and reasonable model for the screened Coulomb 
repulsion V(r, r s ) in the uniform electron gas: he took the sphere of volume n^ 1 as 
the boundary within which the screening charge density is ne and outside of which 
it is zero. We thus have the electrostatic potential due to a point charge — e in the 
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origin plus a sphere around it of radius r s and of uniform positive charge density 
ne. In Hartree atomic units: 

V(r,r s )= i + s) 
V(r,r,)= 

This is equivalent to assuming that in the interacting gas the probability of finding 
three electrons in a sphere of radius r s is exactly zero, an assumption which is nearly 
true. In fact, numerical estimates of this probability for the electron gas show that 
it is indeed small. 15 (At r s = 5 the ratio between the probabilities of finding three 
and two electrons in the same sphere of radius r s is about 1/11; for larger r s this 
ratio is lower, and for smaller r s it is higher, being about 1/7 at r s = 0.) Thus, 
for interelectronic distances r < r s we expect the Overhauser potential to be close 
to the true potential felt by an electron moving in a uniform electron gas when 
another electron is fixed at the origin. In the region r > r s the potential is set 
to zero, and so is not expected to be reliable. We also expect to have results that 
become more accurate as the density decreases, since the probability of having 
three electrons in the same sphere of radius r s becomes lower and lower. Finally, 
at high and intermediate densities our results will be much closer to the true g(r) 
for antiparallel-spin correlations than for parallel-spin ones. When two electrons of 
opposite spins are in the same sphere of radius r s , a third electron is excluded from 
the sphere because of both the Pauli principle and the Coulomb repulsion. For a 
pair of parallel-spin electrons, only the Coulomb repulsion prevents a third electron 
of opposite spin from entering the sphere of radius r s , a mechanism which becomes 
less efficient as the density (and thus the kinetic energy) increases. 



r <r s 
r > r s . 



(27) 



5.2. Solution of the model 



The radial Schrodinger equation corresponding to the Overhauser potential 
can be solved analytically (but not in closed- form), 14 and the functions Re(s, q, r s ) 
(where s = r/r s and q = kr s are scaled variables) are found. As an example, in the 
left panel of Fig. 2, the radial functions corresponding to r s = 4, £ = 0, and three 
different values of the scaled relative momentum q = kr s are reported, together with 
the effective potential. For an unpolarized uniform gas, the probability distribution 
for q can be obtained from Eq. (18) by setting hp — ^(^ t ) 1 ^ 3 , 

M = (28) 

with q ranging from to (9-7T/4) 1 / 3 w 1.91916. This p(q) is exact for a noninteracting 
gas. As said, p(q) for the interacting gas slightly deviates from Eq. (28), and depends 
explicitly on r s . However, in the region where the potential is reliable, r/r s < 1, 
we see from Fig. 2 that the g-dependence of Rt is very weak, as already pointed 
out by Overhauser. 13 We thus expect to have no significant change in the short- 
range (r < r s ) part of g if we use an interacting momentum distribution instead of 
Eq. (28). 

In the right panel of Fig. 2 the g(r) obtained by solving the Overhauser model 
are reported, and compared with the newest diffusion Quantum Monte Carlo (QMC) 
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Figure 2: Left panel: the Overhauser potential V(r,r s ) for r s — 4, together with 
the radial wavefunctions Ri = o for different values of the scaled relative momentum 
q = kr s . Right panel: pair-distribution functions g(r) of jellium computed with the 
Overhauser potential compared to the Quantum Monte Carlo data of Ortiz, Harris 
and Ballone. 16 



data of Ortiz, Harris, and Ballone. 16 In the range < r/r s < 2, a truncation of 
the infinite sum over t in Eqs. (25) and (26) at £ max = 7 is enough to reach good 
convergence. 14 We see that there is accurate agreement with the QMC data for 
0.5 < r/r s < 1 for a wide range of electron densities. In the shortest-range region, 
r/r s < 0.5, the QMC data are known to suffer large errors and are not so reliable. 
(In this region there is in fact a significant discrepancy between the data from Ref. 16 
and those from Ref. 17.) Thus, for r/r s < 0.5 the present treatment should provide 
results much closer to the true g(r). As said, for r/r s > 1 the results obtained with 
the Overhauser potential are not reliable. 



6. SUMMARY AND CONCLUSIONS 

Two-electron wavefunctions (geminals) for the uniform electron gas of density 
n = 3/47rr;? are found using the simple screened Coulomb potential proposed by 
Overhauser. These wavefunctions give pair-densities in agreement with the QMC 
simulations in the short-range (r < r s ) region. They should thus be accurate for 
r < r s , and could be used as a comparison for other two-electron approaches to 
the jellium model, and for testing reduced-density matrix energy functionals. In 
particular, it would be of great interest to compare them with the geminals which 
will be computed following the approach described by P. Zicsche in this book. 

The main lack of the present treatment is the violation of the particle-conservation 
sum rule. A possible solution to this problem is to write down a self-consistent set of 
equations in which the effective potential is unknwon, but the "exact" pair-density 
of jellium 18 ' 19 is used to generate it. 
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